Electron Evaporation from an Ultracold Plasma in a Uniform Electric Field 
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Electrons in an expanding ultracold plasma are expected to be in quasi-equilibrium, since the 
collision times are short compared to the plasma lifetime, yet we observe electrons evaporating out 
as the ion density decreases during expansion. We observe that a small electric field that shifts the 
electron cloud with respect to the ions increases the evaporation rate. We have calculated the spatial 
distribution of a zero-temperature electron cloud as a function of applied field and ion density, which 
is assumed to be Gaussian at all times. This calculation allows us to predict the flux of cold electrons 
from the plasma at all times, and is in good agreement with our observed electron signal. 

PACS numbers: 32.80.Fb, 52.50.Jm, 52.55.Dy 



Ultracold neutral plasmas (UNPs) [T|, Q are created by 
the photoionization of laser-cooled neutral atoms. They 
have typical densities of 1 x lO^cm^-^ and electron tem- 
peratures ranging from 1-lOOOK. UNPs are typically un- 
confined and expand due to electron pressure into the 
surrounding vacuum with lifetimes of a fevi^ hundred mi- 
croseconds, limited by stray electric fields. During this 
expansion, the electron temperature is changing dynam- 
ically from the competition of adiabatic and evaporative 
cooling with various heating processes, such as three- 
body recombination. Previous experimental work has 
focused on measuring the global electron temperature 
[^-Q, but the thermal distribution of the electrons has 
also drawn interest ^-(9(1 . Since the plasma electrons 
are held in a finite-depth potential well, a true Maxwell- 
Boltzmann distribution cannot exist, leading to predic- 
tions of a Michie-King distribution Q , which has been 
developed for gravitationally bound systems in astro- 
physics. This is in contrast to a truncated Boltzmann 
distribution that has been proposed for a similar sys- 
tem in quasi-equilibrium, trapped ultracold atoms un- 
dergoing forced evaporation for the production of Bose- 
Einstein condensates [l3|- Electrons in UNPs may also 
provide insight into other physical systems, as their be- 
havior shares many similarities with electrons moving in 
laser- irradiated atomic clusters [11] and stars in globular 
clusters 

Probing electron dynamics is difficult in UNPs since 
the low total number of electrons prevents the use of 
electric probes inside the plasma. Instead, electrons are 
detected by applying a small uniform electric field across 
the plasma to direct escaping electrons to a charged par- 
ticle detector. Electron detection has been used to mea- 
sure collective oscillations of the plasma [H, [T 3| . Ryd- 
berg atom formation 
electron temperature 



|, plasma instabilities [15j, and 
All of these measurements 
involved applying perturbing fields, from dc pulses to mi- 
crowaves, and observing the change in the electron emis- 
sion. Electron detection has also been a main diagnostic 
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in studies of the spontaneous evolution of Rydberg gases 
into plasmas [l^ [l3| and of a UNP in a strong magnetic 
trap [18| . But an explanation of the main feature of the 
unperturbed electron emission signal is missing from the 
literature. 

A typical emission curve from our experiment over the 
entire lifetime of the plasma is shown in Fig. [TJa) , show- 
ing a prompt peak of electrons emitted during plasma for- 
mation, followed by electrons continuously emitted dur- 
ing the expansion phase. In principle, this signal contains 
information about the electron distribution and evapora- 
tion processes. At a given time, each electron is bound 
to the plasma in a potential well of finite depth, created 
by the combination of the slow-moving ion cloud, the 
applied electric field and the other electrons. Since the 
electron collision times are short compared to the plasma 
expansion, electrons will evaporate out of the plasma un- 
til a quasi-equilibrium is established. This evaporation 
differs from that for neutral atoms in that the trapping 
potential depends on the electron number; for each elec- 
tron lost, the plasma potential well binding the remain- 
ing electrons deepens. As the density decreases from ex- 
pansion, the external field can penetrate farther into the 
plasma, spilling more electrons and forcing the evapora- 
tion to continue. 

In this paper, we separate the effects of the external 
electric field on electron emission and find that this con- 
sideration alone can reproduce the basic shape of the 
electron emission curve. We use a zero temperature ap- 
proximation, and present an electrostatic calculation that 
gives the electron spatial distribution for a fixed ion den- 
sity and applied field. This gives the maximum number 
of electrons that can be held in the plasma in the ab- 
sence of evaporation. Applying this calculation at all 
times during expansion allows us to predict the flux of 
cold electrons from the plasma, which we flnd is in excel- 
lent agreement with the observed electron signal. 

The creation of an ultracold plasma in our experiment 
has been described in previous work fl^. Briefly, we 
prepare a sample of 5 x 10^ metastable Xe atoms in 
a magneto-optical trap. The atomic density is roughly 
Gaussian with an rms radius of « 300/irn. We create a 
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FIG. 1. (a) An example of our canonical electron current 
signal, (b)-(d) The integrated data signal (black solid) com- 
pared to theory (blue dashed) with Eext ~ 50 mV/cm. Neu- 
trality curve in (b) corresponds to the electron signal in (a). 
Plasma parameters are (b) Ni — 1.1 x 10'',i;o ~ 47m/s, (c,i) 
Ni = 0.23xl0^•^;o = 75m/s, (c,ii) Ni = l.OxlO^.wo = 70m/s, 
and (d) Ni = 0.92 x lO", -uo = 95m/s. Electron signals are the 
average of 50 shots with characteristic statistical uncertainty 
shown in (b). 



plasma using two-photon photoionization with one pho- 
ton at the cooling transition (882 nm) and the second 
from a tunable pulsed dye laser (514 nm). The initial 
energy given to the electrons, E^., is controlled by tun- 
ing the energy of the pulsed laser above the ionization 
limit. After creation, the plasma loses a few percent of 
the electrons until a large enough charge imbalance exists 
to trap the remaining electrons, which quickly thermalize 
(the electron thermalization time is w 20 ns at the cen- 
ter of the plasma at formation densities) and establish 
a quasi-equilibrium. This initial loss is represented by 
the prompt peak in the signal in Fig. [TJa). During the 
expansion phase, the ion distribution maintains a Gaus- 
sian density following a self-similar expansion described 
by cr2(i) = 0-2(0) + y^t'^ with a rate of vo = 50- lOOm/s, 
driven by the thermal electron pressure. 

Two wire mesh grids are located 1.5 cm on either 
side of the plasma. An electric field (typically just 5- 
10 mV/cm) apphed between these grids directs electrons 
out of the plasma region, through a third grid and onto a 
microchannel plate detector. By controlling the voltages 
on the grids and detector, we can detect either electron 
or ion emission from the plasma. A schematic showing a 
plasma polarized by this field is shown in Fig. (^fa). 



As the plasma expands, the electron and ion densities 
decrease and the electron temperature decreases follow- 
ing adiabatic expansion. Previous measurements show 
the electron temperature scales approximately as t~'^^^ 
for times > 15^s Q. The combined effect means the 
thermalization times for electrons steadily increase but 
remain less than a few /is over the entire plasma lifetime. 
This allows us to assume that the electrons maintain a 
quasi-equilibrium within the continually decreasing po- 
tential well. 

For our calculation, we first consider a fixed point in 
time and allow the electrons to move and establish equi- 
librium within the stationary ion density. An example of 
a Gaussian ion cloud tipped by a uniform external field is 
shown in Fig.^b) (dashed curve). We treat the electrons 
as a zero-temperature fluid, and seek the maximum num- 
ber of electrons such a potential can hold and the result- 
ing spatial distribution. The zero-temperature approxi- 
mation requires a net zero electric force on all electrons 
in equilibrium, so the maximum number of electrons can 
be determined by adding electrons until the combined 
electron and ion potential becomes flat (Fig. [UJb) (solid 
curve)). In the absence of an external field, the solution 
is an electron density that exactly matches the ion den- 
sity everywhere in space. But an external field reduces 
the number of electrons that can be held, pushing the 
extracted electrons away from the system. The electron 
density will still exactly match the ion density in some 
region of space with a sharp boundary and be zero out- 
side that boundary. The shape of the boundary is defined 
such that the ion density outside the boundary has the 
correct shape to exactly cancel the external field every- 
where within the electron cloud, resulting in a net zero 
force on all electrons. 

To find the proper shape of the electron cloud, we use 
a numerical algorithm first used by B. Breizmand and A. 




FIG. 2. (a) A schematic of our setup including a contour 
plot of ion density and electric grids. The solid line shows 
the sharp edge of the cold electron cloud, and the dashed 
line is the 3cr radius of ion cloud, (b) The potential along 
the symmetry axis of the ion distribution plus external field 
(dashed) and the total distribution including cold electrons 
(solid). 



3 




FIG. 3. A potential map of the plasma showing lines of 
equipotential for a — 0.0035, which corresponds to an electron 
fraction of 0.47. The thick solid lines are the boundaries of the 
ion and electron clouds. The potential is constant everywhere 
inside the electron cloud boundary. 



Arefiev to analyze a similar problem in laser-irradiated 
cluster physics [2^, |2lj. We define a boundary to the 
electron cloud that is initially at the edge of the ion dis- 
tribution. We take the ion density to be a Gaussian with 
a cutoff at 3 cr. The cutoff reflects the predicted existence 
of an ion density spike near the 3 a point of the expand- 
ing plasma that arises from the initial charge imbalance 
d, . We calculate the electric field at the boundary 
and then displace the boundary by an amount propor- 
tional to the local value of the field. After repeating this 
process over many iterations, the boundary motion slows 
down and eventually stops as it reaches the case of zero 
field at all boundary points. Inside the boundary the 
electron density exactly matches that of the Gaussian 
ions and the integration of this gives the total electron 
number. An example of the electron boundary shape is 
given in FigO 

We perform this calculation at all times during expan- 
sion, assuming the ion cloud center is fixed and the size 
follows the ideal self-similar expansion. As the density 
decreases, the electric field has a stronger polarizing ef- 
fect such that fewer and fewer cold electrons can be held 
in by the ions. Eventually the density is low enough that 
the electric field dumps all electrons. We calculate the 
maximum number of cold electrons as a function of time 
(ion density). 

The shape of the boundary and the fraction of electrons 
to ions, Ne/Ni, depends only on the ratio of the external 
field, Eext, to the characteristic field of the ions. We find 



Ne/Ni for one set of dimensionless values of 

where eo is the electric constant, e is the electron charge, 
and Ni is the number of ions. Measurements of the 
plasma expansion velocity and density are sufficient to 
compare to experiment. 

Figure [TJb) shows the fraction of electrons in the 
plasma as a function of time, found by integrating the 
signal in Fig.[TJa), and the corresponding calculated cold 
electron fraction. The calculation predicts the general 
shape and timescale of the signal over a wide range of 
initial parameters. As expected, the calculation overes- 
timates the plasma neutrality during the prompt loss of 
thermal electrons, the discrepancy greater for higher Ee 
and lower Ni [l^ . The expansion velocity of the plasma 
is measured using the plasma response to an rf field , 
and the number of plasma ions is left as a fit parameter. 
Fitted Ni values fall within the uncertainty of previous 
independent measurements made using MOT diagnostics 
and other plasma processes. 

Figure m shows a direct comparison of data taken with 
different strengths of electric field. Stronger fields shorten 
the plasma lifetime, but plotting the signals using our 
normalized coordinate a shows a similar rate of elec- 
tron loss for all. Electron signals taken for varied ini- 
tial density and initial electron energy are also found 
to follow this universal curve. For initial energies below 
Ee/kb ^ 60K the shape of the observed signals does not 
match as well to the theory. At these low temperatures, 
an increased three-body recombination rate converts as 
much as 15% of the plasma ions into Rydberg atoms, 
@j (our calculation has assumed a constant ion num- 
ber) and the model of self-similar Gaussian expansion is 
less accurate [4]. 

Late in its lifetime, the plasma is large enough for the 
edges of the ion cloud to start crossing the mesh grids. 
By monitoring the ion current signal, we determine that 
the first ions start crossing the grid when the grid po- 
sition is between 3 a{t) and 4 cr(t). The ion signal also 
shows a sharp turn-on, indicating a sharp edge to the 
ion density, which supports our assumption of a Gaus- 
sian density with a sharp drop at «3 a. For times after 
this, a more correct density includes a planar cut-off at 
the locations of the grids. In Fig.HJb), we compare only 
the times before the ions start crossing the grid. Prac- 
tically this makes little difference, as the grids only cut 
off the low density outer regions of the plasma. We have 
performed the cold electron calculation with the grids in- 
cluded for a large and slow (9 x 10^ ions and t^o — 60 m/s) 
plasma with a low electric field (15 mV/cm), where the 
lifetime is long. The result gives lower neutrality for all 
times after the plasma starts crossing the grid, but the 
differences are less than 5 percent compared to the case 
where we ignore the grids. 

Above 200 mV/cm, the signal scaled with a begins to 
diverge from the universal curve seen in Fig. Hl^b). For 
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FIG. 4. (a) Integrated electron signals for varied external 
field, and (b) the same signals plotted vs the normalized a. 
The black dashed line is the theoretical result. Data signals 
are the average of 50 shots with -Be/fci, = lOOK and A'^; = 
1 X 10®. The ion cloud starts crossing the grids at 71/is. 



strong fields, the neutrality drops quickly enough that it 
alters the ion expansion, both from the linear acceleration 
of free ions in the applied field and from the Coulomb re- 
pulsion between the ions. The Gaussian density approx- 
imation then becomes invalid. Residual electric fields 
in our vacuum chamber prevent a good fit down to ar- 
bitrarily small field, and ultimately limit the maximum 
observable plasma lifetime. Stray fields along the grid 
axis can be canceled with slight adjustments to the grid 
voltages. We monitor the total number of electrons de- 
tected over the full lifetime and see a drop from all elec- 
trons detected to none as one grid voltage is changed over 
a 25 mV range. The midpoint of this transition should 
be close to zero axial field and typically comes with a 
20-30 mV difference between the grids. This difference 



is taken into account in our reported field values. Stray 
transverse fields are not canceled, but we estimate them 
to be no greater than 20 mV/cm based on the quality of 
our fits near zero axial field. 

Expanding ultracold neutral plasmas are systems in 
quasi-equilibrium, sharing this aspect with a wide vari- 
ety of physical systems from atomic to galactic clusters. 
We have presented data of electrons evaporating from the 
expanding plasma that is in excellent agreement with a 
zero-temperature electron model, where the electrons ar- 
range themselves in an asymmetric density distribution 
to exactly cancel the ion plus external electric fields. This 
result does not directly address the question of the form 
of the electron energy distribution function, but rather 
serves as a foundation that can be used in further exper- 
iments. Experiments using the electron leakage to probe 
electron energies must find a way to separate or eliminate 
the spilling of electrons from thermal evaporation. 

Eliminating the effect of the external field might be 
done by dumping a significant percentage of the electrons 
at early time to create a deep potential well, and watch- 
ing for electrons that escape in the following time interval 
when the cold electron flux should be zero. Alternatively, 
we can imagine decreasing the electric field continuously 
to counter the effect of expansion and prevent any cold 
electron loss (i.e. keep a constant) and detect the emit- 
ted thermal electrons. Our ability to directly detect the 
emitted particles as well as perturb the potential in a time 
dependent manner suggests that UNPs are a good can- 
didate system to address fundamental questions about 
systems in quasi-equilibrium. 
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